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Abstract. Infrared emission from radiatively heated dust possesses general 
scaling properties. The spectral shape is independent of overall luminosity 
when the inner boundary of the dusty region is controlled by dust subli- 
mation; the only relevant property of the heating radiation is its frequency 
profile. For a given type of dust grains, the emission from red giants and 
supergiants is essentially controlled by a single parameter — the overall 
dust optical depth. This leads to tight correlations among different spectral 
properties that explain many available observations and enable systematic 
studies of large data bases. 

1. Introduction 

The radiation received on earth from late-type stars has undergone signif- 
icant processing in the surrounding dust shell. Interpretation of the obser- 
vations therefore necessitates considerable theoretical effort, involving de- 
tailed radiative transfer calculations. These calculations have traditionally 
required a large number of input parameters that fall into three categories: 

The Star. Stellar input properties include the mass M*, luminosity L*, 
temperature T* and the mass- loss rate M. 

The Shell. This is described by a density profile p(r), defined between 
some inner and outer radii r\ and r<i, respectively. 

Dust Properties. The dust abundance can be expressed via the dust-to- 
gas ratio pd/p- The properties of individual grains must be specified, too, 
and quantities widely employed include the grain size a, the solid density p s , 
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sublimation temperature T su b and the absorption and scattering efficiencies 
Q abs and Q sea- 

Once these input properties are prescribed, they are plugged into a de- 
tailed radiative transfer calculation whose output can be compared with the 
observations. The most widely used output quantity is the spectral shape 



where F\ and F are the observed flux density and bolometric flux, respec- 
tively. Another quantity that can be compared with observations when the 
angular resolution is sufficiently high is the surface brightness. 

The rather large number of input parameters that had to be specified 
in traditional calculations creates two major practical problems. First, the 
volume of parameter space that must be searched to fit a given set of 
observations can become prohibitively large. Second, and more serious, even 
when a successful fit is accomplished, its uniqueness is questionable and the 
model parameters cannot be trusted as a reliable indication of the source 
actual properties. 

2. Scaling 

Much of the input required in past calculations is in fact redundant. In 
Ivezic & Elitzur 1995 (IE95 hereafter) we were able to show that the IR 
emission from late-type obeys general scaling properties that greatly reduce 
the number of input parameters required to fully specify the radiative trans- 
fer problem. The primary input involves only the dust, for which we must 
specify two types of properties: (1) Optical properties, specified through the 
normalized spectral shape of the extinction efficiency q\ = Q\/Q\ , where 
Ao is some arbitrary fiducial wavelength; this spectral shape is controlled 
by the grain chemistry and size. (2) The overall optical depth r^; this is 
controlled primarily by the mass-loss rate. Additional input properties have 
only secondary significance and involve the dust sublimation temperature, 
Tsui,, the relative thickness of the dust shell ti/tx and the stellar tempera- 
ture Note that T* is the only stellar property that enters (and only in 
a secondary role). In particular, the stellar luminosity is irrelevant. 

The proof of scaling is quite simple. The shell inner radius r\ is con- 
trolled by dust sublimation, namely, T^(r{) = T su b. Introduce the dimen- 
sionless radial distance y = r/r±. Then the radiative transfer equation 
becomes 




(1) 
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for y > 1. In this form, the equation contains no reference to either densities 
or dimensions. The only scale is determined by r T , geometrical quantities 
enter only through the dimensionless radius y. The density enters only 
through rj = n(y)/J n(y)dy, its dimensionless, normalized radial profile. 
Because the outflows around late-type stars are controlled by radiation 
pressure on dust grains, the profile rj is uniquely determined by t t when 
the effects of gravity and dust drift are negligible (IE95). In fact, under 
these circumstances, the analytic expression 



provides an excellent approximation to the actual density profiles we find 
in our detailed numerical calculations. This function requires as additional 
input the ratio of initial to final velocity, fi/foo- However, this parameter 
has a negligible effect on the result for its typical values, < 0.1. 

Additional input required is the stellar radiation. This can be specified 
by its flux density F*\ = x f*x, where is the bolometric flux. Each 
value of Fi = F*(y=l) uniquely determines a corresponding value of T\, 
the dust temperature at y = 1, so the reverse is also true. And since the 
inner boundary is controlled by dust sublimation, T su f> fixes the value of T\ 
and through it of F\. When the stellar luminosity varies, the shell adjusts 
its inner boundary so that F\ and T\ (= T su b) stay the same. The radiative 
transfer equation is oblivious to these changes because the inner boundary 
always corresponds to y = 1. 

These scaling properties can be generalized to arbitrary geometries and 
density distributions (Ivezic & Elitzur 1996). The spectral shape of dust IR 
emission is independent of overall luminosity when the inner boundary of 
the dusty region is controlled by dust sublimation; the only relevant prop- 
erty of the heating radiation is its spectral shape. Densities and geometrical 
dimensions are likewise irrelevant; they enter only through one independent 
parameter, the overall optical depth. The geometry enters only through an- 
gles and aspect ratios. Dust properties enter only through dimensionless, 
normalized distributions that describe the spatial variation of density, and 
the wavelength dependence of scattering and absorption efficiencies. 

3. CONSEQUENCES 

Scaling implies that for a given type of dust, the spectral energy distribu- 
tions (SEDs) of late-type stars are controlled almost exclusively by overall 
optical depth and thus form a one-parameter family. Therefore, a single 
point on the normalized spectral shape f\ determines the entire function 
and any two spectral properties should be correlated with each other. Fig- 
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Figure 1. Test of correlations predicted by scaling. Open squares mark the data points 
for O-stars, solid circles C-stars. Curves describe the model predictions for various grain 
compositions, as marked. 



ure 1 displays in its left panel the distribution of normalized fluxes at 12 
and 60 /im for a sample of 89 stars. These are all the IRAS objects identified 
as late-type stars for which we were able to find listings in the literature 
for both total fluxes, luminosities, terminal velocities and mass-loss rates. 
As expected from scaling, the normalized fluxes do display a tight corre- 
lation that holds over the entire observed range, covering more than three 
orders of magnitude. The spread can be attributed to the differences in dust 
properties between carbon- and oxygen-rich stars. The right panel displays 
the corresponding distribution of {vF p )qq/F and [25]— [12] color. The ex- 
pected correlation again is evident, the separation between C and O stars 
is more pronounced. The lines are our computed model spectra, displaying 
a close agreement with the data for various relevant chemical composi- 
tions (Sil stand for astronomical silicate; am.C for amorphous carbon; Oli 
for olivine). The agreement between model predictions and observations is 
better than factor 2 in almost all cases. In addition, O-stars and C-stars 
clearly separate and congregate in accordance with the trend of the model 
curves for corresponding chemical compositions. This behavior is particu- 
larly prominent in the right panel. 
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Figure 2. IRAS color-color diagrams of all sources in the VH-window with reliable data. 
Solid lines represent theoretical sequences for different grain compositions, as marked. In 
all mixtures, the abundances ratio of the first to second component is 1:4. In the lower 
panel, the order of the tracks from bottom to top around [25]-[12] ~ is graphite, Oli. 
+ Sil, Sil, Sil + am.C, am.C and SiC + am.C. 



Location in IRAS color-color diagrams has become a widely used indi- 
cator of the nature of a source. Van der Veen & Habing (1988) identified 
the appropriate IRAS region for late-type stars as [60]-[25] < and [25]- 
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[12] < 0.6, which we will refer to as the VH-window. Figure 2 displays the 
color-color diagrams for all IRAS sources with reliable data whose colors 
fall in this window. Our model predictions for selected grain compositions 
are displayed as tracks. Position along each track is determined by t t and 
all tracks originate from the same spot, r T = 0, corresponding to a black- 
body spectrum at 2500 K convoluted with the IRAS instrumental profiles. 
Distance from this common origin along each track increases with t t . As r T 
increases, the dust emission becomes more prominent and the tracks of dif- 
ferent grains branch out. The tracks for purely Si- and C-based grains out- 
line the distribution boundaries of most IRAS sources in the VH-window, 
verifying the VH conjecture that this is the color-color location of late-type 
stars. The data points fill the entire region between these tracks and if this 
spread is real, it requires a similar spread in the optical properties of the 
grains. This could reflect chemical mixtures, as we proposed in IE95, or 
perhaps grain impurities. 

It is important to note that we have removed from our sample sources 
whose IRAS fluxes are severely contaminated by cirrus emission. This con- 
tamination is measured by the cirrus contamination index 

where cirr3 is the 100 //m cirrus flux at the location of the source and F(60) 
is its listed 60 fim flux. In IE95 we show that the [100]-[60] color and CCI 
are perfectly correlated in sources with CCI > 2, indicating that the IRAS 
60 and 100 fim fluxes of these sources are unreliable. Sources in the upper- 
left corner of the color-color diagrams were thought to provide evidence 
for cool shells (Willems & de Jong 1986), but almost all of these sources 
turn out to be cirrus contaminated. Of the total 292 sources displayed in 
figure 2, only 11 late-type stars (4%; marked with open squares) can be 
identified as a group whose colors are not well explained by steady-state 
radiatively-driven winds. Furthermore, among these 11 sources, 10 have 
borderline contamination with 1 < CCI < 2. Based on new simultaneous 
observations in a wide spectral range, Miroshnichenko et al. present in a 
poster paper model fits for 4 carbon stars. The results are reproduced in 
figure 3, showing an excellent agreement with the model for steady-state 
winds without any additional shells. The only parameter that varies from 
model to model is the visual optical depth t v . Two of the displayed stars, 
RY Dra and VY UMa, belong to the group marked with squares in figure 2. 
They are borderline contaminated and show a 100 pLvn excess proportional 
to CCI, so clearly they do not give any indication of detached shells. There 
is only one clearly uncontaminated C-star, U Ant, in the "cool shell" zone. 
Indeed, this star has a thin CO shell (Olofsson et al. 1990). We conclude 
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that, because of potential cirrus contamination, the SED is unfortunately 
not a good indicator of detached shells. 
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Figure 3. Model fits for four carbon stars with new ground based observations (triangles; 
Miroshnichenko & Kuratov 1996). Circles denote LRS data, squares are IRAS PSC and 
RAFGL data. The only adjusted model parameter is the visual optical depth t v . 



Thanks to scaling, the SEDs of large samples can be modeled system- 
atically. We have embarked on a large-scale modeling project in which we 
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are fitting all IR data of all the late- type stars in the IRAS catalogue (~ 
5,000 sources). The outcome of this project will be a supplemental cata- 
logue, listing the optical depths and dust properties of the largest sample 
yet modeled with a single, consistent approach. 

4. CONCLUSIONS 

The structure displayed in infrared color-color diagrams is a result of general 
scaling properties of radiatively heated dust. Sources segregate into families 
that differ from each other because of dust properties, not because of the 
central source. Within family, location of each source is controlled by its 
optical depth. Concerning late-type stars — properly selected IRAS data is 
almost fully explained by steady-state winds. The spectral shapes of sources 
free of cirrus contamination do not give evidence for detached shells in more 
than a handful of sources. 
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gratefully acknowledged. 
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